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We derive the leading asymptotic limit of the Wigner 3j-symbol from a stationary 
phase approximation of a twelve dimensional integral, obtained from an inner product 
between two exact Bargmann wavefunctions. We show that, by the construction 
of the Bargmann inner product, the stationary phase conditions have a geometric 
description in terms of the Hopf fibration of into R'^ x R-^ x M^. In addition, 
we find that, except for the usual modification of the quantum numbers by 1/2, 
the imaginary part of the logarithm of a Bargmann wavefunction, evaluated at the 
stationary points, is equal to the asymptotic phase of the 3j-symbol. 
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I. INTRODUCTION 



This paper presents a new derivation of the asymptotic hmit of the Wigner 3j-symbol, 
using stationary phase approximation on an inner product of wavefunctions in the Bargmann 
representational"-. The asymptotic hmit of the 3j-symbol was first derived in 1968 by Pon- 
zano and ReggeAl. Since then, there were several subsequent derivationsi'^'^ii^i^'i^. Despite 
the long history behind this result, this paper uncovers some new connections between the 
asymptotic formula and the Bargmann representational"-, which is an exact representation 
of SU (2) by holomorphic functions. 

We find that, by construction of the inner product in the Bargmann representation, the 
stationary phase conditions are formulated in terms of the projections of three vectors in a 
triangle through the Hopf fibration. In addition, the imaginary part of the logarithms of the 
Bargmann wavefunctions is equal to the asymptotic phase of the 3j-symbol, albeit, without 
the usual modification of the quantum numbers by 1/2. Finally, we point out that the asymp- 
totic wavefunctions appearing in the recent study of the Wigner 3nj-sjmho\a^^'—^^^'^^^ are 
actually Bargmann wavefunctions in disguise. 

We now give an outline of the paper. In section [TTl we give some background on the 
Bargmann representation, and write down the exact Bargmann wavefunctions for the eigen- 
states in the definition of the 3j-symbol. In section llllt we express the 3j-symbol as a 
twelve- dimensional integral. We then apply a stationary phase approximation to this inte- 
gral to derive the asymptotic limit of the 3j-symbol in section HVl After the derivation of the 
asymptotic formula, we relate the Bargmann wavefunctions to the holomorphic wavefunc- 
tions used in geometric quantization in section |Vl and to and the multidimensional WKB 
wavefunctions in section IVII The last section contains some comments and conclusions. 

II. THE BARGMANN WAVEFUNCTIONS 

Bargmann constructed a unitary representation of the quantum rotation group SU (2) in 
a series of papers^"- in the 1960s. Bargmann's construction used Schwinger's construction 
for the generators of SU{2) from two sets of commuting boson operators. Instead of using 
the quantum harmonic oscillators, however, Bargmann implemented the boson operators in 
terms of multiplication and differentiation in a complex coordinate. The result was a Hilbert 
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space that consisted entirely of holomorphic functions. 

We now briefly describe the construction of the Bargmann representation. Let Zf^, /i — 
1, 2, denote coordinates on C^. The bosonic commutation relations, 



[...^1 = 1. (1) 

allow us to form two independent sets of boson creation operators z^ and annihilation 
operators , jj, = 1,2. Each set of boson operators has a number operator. 



= z,d,^ , (2) 

whose eigenvalues are integers, and whose eigenfunctions are monomials in z^. Following 

Schwinger's construction of the generators of SU{2) in terms of boson creation and annihi- 
lation operators, we define the SU (2) generators in the Bargmann representation by 



Ji = ^M^Mi^T^ ' (3) 



d_ 

' dzu 

where o"* are the Pauli matrices. These generators satisfy the usual SU{2) Lie algebra 
commutation relations. 



[Jj, Jj\ = ieijkJk ■ (4) 



The Casimir operator is 



= 7(7 + 1), (5) 

where 7 is given by 

7=i(7Vi + 7V2). (6) 

Since the eigenvalues of iVi, N2 are integers, the eigenvalues of 7 are half- integers j = (ni -|- 
n2)/2, and the eigenvalues of are given by j{j + 1). 

The Hilbert space for a single Bargmann space, T2, consists of holomorphic functions on 
C^. That is. 



•F2 = {/eC~(C2)| (/|/)<oo}, (7) 



where the inner product is given by 



{f\g)= / f{z)g{z)'^d'z, (8) 

and where d"'z denotes the usual Euchdean measure on = R^. 

In the recouphng theory of three angular momenta, the Hilbert space is the tensor product 
of three copies of the Hilbert space for a single Bargmann space. The wavefunctions are 
holomorphic functions on = x x C^. Let us denote the coordinates by (-Zr/^) = 
(-211, -212, -221, -222, -231, -232), r = 1, 2, 3, /i = 1, 2. Let us denote the operators in Eq. ([3]) and (E]) 
that act on the rth angular momentum space by Ir and Jri, respectively. 

We now find the eigenstates that appear in the definition of the 3j-symbol. The basis state 
i^jmiz) = {z\iimij2m2i2,'m-i) satisfies the following eigenvalue equations and normalization 
condition, 



Ir1pjm{z) = ^i^rl + Nr2)^jm{,z) = jri'jmi.z) , (9) 

Jrzi^jmiz) = ^i^rl " Nr2)lpjm{z) = mrtpj„,{z) , (10) 
ii^jml^jm) = 1 , (11) 

for r = 1, 2, 3. The eigenvalue equations determine the degrees of the monomial in Zrfj,, and 
the normalization condition fixes the constant in front. We find 

^a^u Z2) = n , ^ z^"^'- ^r"^^ ■ (12) 

We now find the rotationally invariant wavefunction ipmviz) = (-2IJ1J2J3O), which satisfies 
the following eigenvalue equations and normalization condition. 



Iri^in^iz) = ^{Nrl + Nr2)i'mv{z) = >z/'inv(^) , (13) 

Jilpin^iz) = Jr}j ^/'inv(-2) = , (14) 

(V'invl^inv) = 1 , (15) 



for r = 1, 2, 3, and i = 1, 2, 3. 
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We use the diagonal SU{2) action generated by the generators Jj to find ipirwiz). Let 
us denote the diagonal group action by T{/, where U G SU{2). Its action on a Bargmann 
wavefunction f{z) is given by 

Tu{f{z)) = fiUl^z,,, Ulz^,, Ulz,,) . (16) 
Since U is unitary, the three determinants 

^1 = Z2lZ^2 — ^31^22 ; ^2 = ^31^12 ~ ZiiZ^2 ; ^3 = 2;il2;22 ~ 2;2l2;i2 , (17) 

are invariant under Tu. Using these determinants to construct ipim/{z), we find 

^^^^(^) = , (18) 

where the non- negative integers ^1,^2,^3, which are chosen to satisfy the eigenvalue equation 
in Eq. (fT3i) . They are given by 

ki = 32 + k - ji , h = is + ji - 32 , h = 3i + 32 - 33 ■ (19) 

III. AN INTEGRAL REPRESENTATION OF THE 3j-SYMBOL 

Taking the scalar product between ipjm{z) from Eq. (fT2l) and V'mvl-z) from Eq. (ITSll . we 
obtain an exact integral representation of the 3j-symbol, 



3i 32 33 
nil m2 



{3imi32m2hm^ Ijihh^) (20) 

AT I rl^2 .y^h+rni —ji-mi —j2+m2 —j2-m2 —ja+rn-i —j3-rn3 -|2p 
iv I a ^^11 ^12 ^21 ^22 ^31 ^32 ^ 

(^21^32 ~ ^31^22)*^^ (2:31^12 ~ ^32^11)'^^ (^ll2;22 — 2:12221)'^^ ; 



where the constant in front is 



N ' 



7r6[(ji + mi)!(ji - mi)!(j2 + ^2)1(72 - "^2)!(j3 + "^3)!(j3 - m3)!]V2 
1 



[{3i+ 32 + 33V.h\k2\h\Y/^ 
and where /ci, /c2, are given in Eq. ( IT9|) . 



(21) 
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IV. STATIONARY PHASE APPROXIMATION 



We now apply the stationary phase approximation to the integral expression in Eq. ( 120|) 
in the limit that ji for i = 1, 2, 3 are large. The basic formula for stationary phase approxi- 
mation is given by 

^ V-Hessp(/) 

where p stands for the stationary phase points that satisfy the twelve stationary phase 
conditions 

d.,.Jiz,z) = 0, d,^J{z,z) = 0. (23) 

In the denominator of Eq. fl221) . HesSp(/) is the determinant of second derivatives of / 
evaluated at p. In the general case that the stationary phase points p are not isolated 
points, the sum in Eq. ( 122|) is replaced by an integral over the set of stationary phase points. 

A. Stationary Phase Points 

In the case of the 3j-symbol in Eq. ( I20l) . the phase function / in Eq. ( I22l) is 

f{z,z) = Imfjj^iz) + \mjjinv{z) - l^r/.^ • (24) 

To find the stationary phase points, we use the fact that the Bargmann representation of 
the generators of SU{2) are also differential operators. That is, the twelve conditions ( l23l) 
on the first derivatives of / imply 

i:f{z,z) = 0, J:JXz,z) = 0, lrf{z,z) = 0, JJXz,z) = 0. (25) 

Here we pick the Bargmann operators that appear in the eigenvalue equations ([9]), ffTOj) . 
( !T3|) . ([H]). This way, we can use the eigenvalue equations to easily find the result of the 
Bargmann operators on the first two terms of /. We rewrite the eigenvalue equations. 
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into the following form, 



i* In 1pj^{z) = jr , J*-, In i)jrn{z) =mr, ir In 1pin^,{z) = jr , Jj In V^inv(2) = . (27) 

Since ipj^(z) is anti-holomorphic, Ir and Jj annihilate Inipj^niz)- Similarly, since tpmv{z) is 
holomorphic, /* and J*^ annihilate Inipi^^lz). 

Finally, we apply these Bargmann operators to the last remaining term Y2rii /■ 
The result is 



— ^r ^ ^ l^t'Atl 2 ^ ^ ^''A' '^^M ' '^'''^ — '^rz ^ ^ 2 ^ ^ ij^ z) ^u^ru 1 (28) 

Ir ^IrY: = I Yl ' J^^J^I: \^r,\' = ^ ^ ^^^^^^ • (29) 
rfj, fi r/i rfiu 

The functions Ir, Jrz, and Jj are functions of the Hopf map Jri : x x — )■ x x M^, 
given by 



Jrx — 2 ^r-^j (<72:)^tj^^ri/ — " (^rl^r2 + ^r2^rl) — R-G (^r-l2;r2) ; (30) 

1 V 1 

•^'■y ~ 2 ^^'^ icry)iiuZru = — {ZrlZr2 — ^r2^rl) = Ini (^rl2;r2) , (31) 

Jrz ~ 77 ^ ^ ^r/^ (c'"2)/iv^ri/ = ^ (l^rl | |^r2| ) • (32) 



2 



Putting together the result of apply the Bargmann operators to /, Eqs. ( 17r|) - ( 12^ . the 
stationary phase conditions become 



-/^ — Jt ; Jtz — TTirp J J'i — . (^33^ 

Using the Hopf map, we can interpret the above conditions in x x M'^, which we 
will call the angular momentum space. These conditions are geometrical conditions on three 
vectors Ji, J2, J3. The conditions Ir = jr state the vectors have lengths jr, r = 1,2,3, 
respectively. The conditions Jj = state that the vectors Ji, J2, J3, put together head to 
tail, form a triangle. The conditions Jrz = iTir fix the z projections of the three vectors. 
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Thus a natural way to find the solutions to Eq. ( 133|1 in is to first find three vectors 
Jr in that satisfy the geometric conditions described above. After the vectors are found, 
we lift each vector up to a spinor in the Hopf fiber above. This procedure is carried 
out in Appendix \^ and is very similar to the calculations in the multidimensional WKB 
approach^, except the length of the vectors in that paper are jr + 1/2 instead of jr- The 
stationary phase point p from Appendix |A] is given by 




e 



^ , _ -*^/2pos/3/2 cos 772/2 - sin sin r/2/2 
^ e~^^/2 gin /3/2 cos 772/2 + e^^/2pQg^/2 sin 772/2 

^ . cos /3/2 cos 771/2 + 6*^/2 sin sin r/i/2 

^ e-*T/2 sin /3/2 cos 7/1/2 - e*T/2 cos /3/2 sin 771/2 

-n/../2^f-/^/2,. (3e) 

sin/3/2 



Here the angles 771, 772, /9, 7 are defined in the Eqs. flAip . flA4l) . and flA6p . 

Once we have the solution p with 7 > 0, we can find another solution p' given by 7 = — 17|. 
These two solutions are not isolated in the solution set. To see that, note any set of three 
vectors related to the projection of p in the angular momentum space by a rotation about 
the z axis will continue to satisfy Eqs. ( l33l) . In addition, multiplication by an overall phase 
for each of the three spinors also preserve the conditions in Eqs. fl33l) . Thus, the stationary 
phase points consist of two disjoint 4-tori, generated by an overall rotation about the 2;-axis 
and three overall phases of the three spinors. 



B. The Asymptotic Formula 

After integrating over the two 4-tori of stationary phase points, and applying the sta- 
tionary phase approximation to the integral for the 3j-symbol along the remaining eight 
directions transversal to the stationary phase points, we find 



Jl 32 J3 

rui m2 7723 



gi Im f{p) gi Im f{p') 

+ 



-Hessp(/) A/HesSp/(/) 



(37) 



Here HesSp(/) is the determinant of an 8 x 8 matrix of second derivatives along the eight 
transversal directions. In Eq. fl37|) . one factor of (27r)^ comes from integrating along the four 
angular directions along the stationary points, another factor of (2tt)'^ comes from doing the 
stationary phase approximation along the eight transversal directions. 

The calculation for the Hessian is straightforward but long, so we leave the details of the 
calculation to Appendix [Bl Inserting the result, Eq. (]B47|) . into Eq. ( 137|) . we find 



- (38) 



mi m2 rris 



2^ ^/A 



z 



^ ' cos 5 + - 



26^/51 V 4 

The function /i has the same functional form as the phase function /, but with jr replaced 
by jr + 1/2. This usual modification of the quantum numbers comes from the Hessian. See 
Eq. (IB471) . Explicitly, 



fi{z, 2) = "^ [{jr + 1/2 + rrir) In Zri + (jr + 1/2 - rrir) \nzr2] (39) 

r=l 

3 

|2 



1=1 r/j. 



In the second equality in Eq. (I38l) . we have defined S = lmfi{p), and used the fact that 
Im/i(p') = — Im/i(p). Here is given in Eq. ( 1B46|) . and is the projected area of the 



triangle onto the xy-plane. Thus, the asymptotic phase 5* of the 3j-symbol is equal to the 
imaginary part of the logarithm of the Bargmann wavef unctions, modulo the modification 
of the quantum numbers jr by 1/2. 
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S = J2^i3r + 1/2 + mr) arg(z,i) + (j, + 1/2 - m,) arg(z,2)] (40) 

r=l 

3 

+ ^(/ci + l/2)arg(5,. 

1=1 

Because Si, i = 1,2,3 are invariant under overall SU{2) rotations, the last three terms can be 
evaluated at any point related to the stationary phase points by an overall SU (2) rotation. 



We choose the reference spinor in Eq. (lA3p . The spinors at this point are all real, so the 
last three terms vanish. This is possible because the three vectors form a triangle and can 
be rotated into the xz plane. The remaining terms evaluated at is given by 



_i / Ji cos /3 — mi cos ?72 \ t -i f cos rji — J2 cos /3 
S = Ji cos : ; + J2 cos ' 



sin?72-^i± / V sinr7iJ2_L 

_i /Jicos/3cos?72 — "^lA _i / Ji cos?72 — mi cos/3 

+-'3 cos ( J . n ■ + mi cos — — — 

Ji sm (j sm ri2 J \ sm p Ji_l 

s ?7i - m2 I 

sin /3 J2_L 

where J^j. = \J J'^ — m^, and Jr = > + 1/2- We now calculate the constant factor 
{2TifNe^^f^py2^ in Eq. ([38]). Here 



_i / J2Cos?7i - m2Cos/3\ , , 

-m2Cos — — , (41) 



U ^ -y -y \kl+l/2 \ |fc2 + l/2 I |fc3 + l/2 

1 2:21 2^32 — ^31^221 1^31^12 ~ ^322;il I | ^11 2:22 ~ ^12^21 1 



We evaluate e^''-^ at the stationary point p. From the conditions \zr\\ = jr + ^r, \zr2\ 
jr — rrir, and mi + m2 + m3 = 0, we find 



Y[ |z^i|^''-+'^'-+l/2|^^^|jV-m,.+l/2 I g-E,,^ M 



-2(il+j2+j3) 



1-^'^ + m^|>+m,+l/2 _ |>-m,+l/2 _ (43) 



The remaining factor, 



1^212:32 — ^31^22!'^^^"'"''^ \Z3iZ12 — Zz2Z\\\'^^^^'^ \Z\\Z22 — Z\2Z2\\^''^'^^^'^ , (44) 

is invariant under the diagonal SU(2) actions, so we can again evaluate it at the spinors in 
Eq. (lASp . The result for the three factors are 
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1^21^32 — 2:31 2:22 



|fel+l/2 



V 4j2 js sin(r7i/2) 



fci+1/2 



|4j2j3 (1 - COS?7i)| 



fci+1/2 



2j2j3(l - 



o 2 2 ^ 2 

Jl J2 J3 - 



2j2j3 



fei+1/2 



^(Jl + J2 + J3)'=^+1/2(J2 + J3 - Jl)^'^ + l/2 



2:31^12 — Z32Z11 



fc2 + l/2 



a/ 4jij3 sin(r72/2) 



k2+l/2 



|4jij3(l- COS 7/2)1'^+^/' 



2jlj3(l - 



il ii il ^ 



2jlj3 



A;2+l/2 



^J{Jl + J2 + j3)^'^+^/'(jl + J3 - J2)'=^+^/^ 



kll^22-^12^2l|'=^+'/' 



V4jij2 (sin(r7i/2) 003(772/2) + sin(?72/2) cos(r7i/2)) 

fc3+l/2 



fc3- 



a/ 4jij2 sm((?7i + r72)/2) 

i/|2jU2(l-cos((7?i+r?2)))| 
7|2jU2(l-cos(r/3))|'=^+^/^ 



fc3+l/2 



2jlj2(l 



J3 ii ^2 ' 

2jlj2 ' 



^ (jl + J2 + j3)'=^+^/2(jl + J2 - J^f^^^^'^ ■ 



Multiplying together the above results for Eq. (143|) and Eq. ( 14^ . we find 



gRe/{p) ^ g-2(ji+j2+J3) 



X ^(j, +^3)ii+i.+i3+3/2fcfi + V2^J. + l/2^.3+l/2 



(jl + - "^i)KJ2 + ^2)!(j2 - »^2)!(j3 + ^3)!(j3 " ^3)! 

(27r)5 

X V (ii + i2 + ky-k^W-hKji + J2 + J3) , 
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where we have used Stirhng's approximation for factorials, n\ ^ (27r)^/^ra"'"'"^/^e Most of 
the factorials cancel those that appear in N in (l2Ti) . We find 



^ ^ (49) 



26 

Substituting Eq. dH]) and f H9|) into Eq. f l38|) . we find the asymptotic formula for the 
Wigner 3j-symbol is 

Jl J2 J3 \ , cos(5 + 7r/4) , 

(50) 



27r|A. 



mi m2 



V. BARGMANN WAVEFUNCTIONS AND GEOMETRIC 
QUANTIZATION 

Recently, Roberts^ and Charles^ provided two independent derivations of the Ponzano- 
Regge formula for the Wigner 6j-symbol using wavefunctions from geometric quantization 
of the group SU{2). We now comment briefly on the connections between the Bargmann 
wavefunctions and the holomorphic wavefunctions constructed from geometric quantization 
of SU{2). The wavefunctions that result from the geometric quantization of the coadjoint 
orbits of SU{2), which are 2-spheres, are holomorphic functions on those coadjoint orbits. 
These holomorphic functions can be obtained from the Bargmann wavefunctions by restrict- 
ing them to a CP^ = S'^ subspace of C^. 

Let us focus on a single angular momentum. We start with the Bargmann wavefunction 
for the basis state \jm), given by 

^im(^l, Z2) = ^ Z^^"^ Z^-"^ . (51) 

V(j -m)\{i + m)\ 
The function ipjmizi, Z2) has the scaling property 

Ipj.nii^Zi, \Z2) = A^^ ^jm{Zl, Z2) , (52) 

for A G C. Thus the Bargmann wavefunctions are completely determined by its values on a 
CP^ = S"^ subspace. Let us choose this subspace to be deflned by + |-22p = 1, and set zi 
real. Let ^ = Z2/Z1 parametrize the coordinate on this subspace. Then the basis functions 
ipjm{zi, Z2) restricted to the section CP^ is given by 
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I / /-\ 1 2 

^^'"^^^ " [(j+m)!(j-m)!]i/2 



]l/2 

(53) 



[(j + m)!(i -m)!]i/2 



[(j+m)!(j-m)!]V2(i + |^|2), ' 

where we have used 1 = + |z2p = + I^H- These basis wavef unctions V'jmlO agree 
with those used in geometric quantization of SU (2) on page 177 of Woodhouse^. 



VI. RELATION TO THE SCHWINGER REPRESENTATION 

Another recent approachi'S to the semiclassical analysis of the Wigner Snj-symbols uses 
the WKB wavefunctions of the harmonic oscillators in Schwinger's model of angular momen- 
tum. It turns out the Bargmann wavefunctions are the exact harmonic oscillator states for 
those WKB wavefunctions in the coherent state representation. To show this, we write the 
eigenstates of the harmonic oscillators in the coherent states basis of the Heisenberg-Weyl 
group"*^*^. The coherent states basis are given by \z) = D{z) |0), where |0) is the ground state 
of the simple harmonic oscillator, and D{z) is an element of the Heisenberg-Weyl group, 
given by 

= exp(— |2;p/2) exp(2;d''") exp(— zd) , (54) 

where z = x + ip, and a and d^ are the usual annihilation and creation operators of the 
simple harmonic oscillator. The coherent state wavefunction 

{x\z) = {ttH)^^^^ exp 

is the displaced ground state of the simple harmonic oscillator centered at {x,p). Using the 
exponential generating function for the Hermite polynomials 

°° z^ 

exp(-;22 ^ 2zx) = XI ~r ' 

we find that the eigenstates \n) of the simple harmonic oscillator written in the coherent 
state basis is given by 
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t){z) = exp 



h 



{xp — px) 



(55) 



{z\n) = 



n' 



Z 



n'\ 




dx{'KK) 



Z' 



(57) 



where we have used the orthonormahty of the eigenstates of the simple harmonic oscillator. 

Schwinger's representation uses two simple harmonic oscillators. The basis states | jim) = 
\nxn'ij is the product of two eigenstates of the two harmonic oscillators, where rii = j + m 
and n2 = j ~ m. Thus, the Schwinger wavefunction for \jm) written in the coherent state 
basis is given by 



which is proportional to the Bargmann wavefunction ipjmiz) = {z\jm). Thus, the Bargmann 
wavefunctions are the exact states of the WKB wavefunctions from the multidimensional 
approach. 



In this paper, we have provided a new derivation for the asymptotic formula of the Wigner 
3j-symbol, using exact wavefunctions in the Bargmann representation. In addition, we have 
pointed out the relationship between the Bargmann wavefunctions and other asymptotic 
wavefunctions used in two other recent approaches to the semiclassical analysis of the Wigner 
Sraj-symbols. 

In the derivation in the this paper, as well as in the other two approaches, the stationary 
phase points have natural geometrical interpretation in terms of classical vectors. The 
Bargmann approach and the geometric quantization approach present easy ways to calculate 
the asymptotic phase, whereas the multidimensional WKB approach provides an elegant 
formula for the amplitude determinants. The Bargmann approach, however, is the only one 
that uses the exact wavefunctions directly. This may make it possible to calculate higher 
order terms in the asymptotic series. 

Appendix A: Finding A Stationary Phase Point 

We now find one stationary phase point that satisfy Eq. (l33ll . Let us start with a standard 
orientation for a triangle formed by Ji, J2, J3 by placing J3 along the 2;-axis, as illustrated 




,i+m j-m 



(58) 



VII. CONCLUSIONS 
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in Fig. [H The angles rjr lie in the range < ?7r < vr. and are complementary to the interior 
angles of the triangle with lengths ji, j2, is- From the law of cosine for triangles, we find 



cos rji 



^1 2 ^2 ^2 ^2 ^2 ^2 



and cyclic permutations. 



J. 



2j2j3 



2jlj3 



J. 




(Al) 



FIG. 1. Visualization of the stationary phase conditions | J^l = jr, r = 1,2, 3, and Ji = 0, i = x, y, z, 
as a triangle with edge lengths ji ■, 32 , js ■ 



In that orientation, the vectors in Euclidean coordinates are given by 



Ji = ji 



^ sin?72 ^ 


y cos J 



^ — sin rji ^ 



32 



\ 





cosr^i 



J3 



J 







(A2) 



We find the spinors on the Hopf fiber above these vectors. Because all three vectors lie 
in the xz plane, we can choose all three spinors to be real. We choose 
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We now apply rotations to the vectors in Eq. f lA2p to ensure their z projections satisfy 
Jtz = ^r, r = 1,2, 3. We do this in two steps. First we rotate the vectors in the x-z plane 
about the ?/-axis by an angle /3, < /3 < n, defined by 

mg = ja cos 13 , (A4) 
so that Jsz = f^s- This rotation is illustrated in Fig. O 

Jz 



J X 









\j2 








1713 



FIG. 2. The vectors after applying the rotation about the y-axis by an angle (3 to the triangle in 
the reference orientation. 



Once we have J32 = ms, we rotate the vectors about the vector J3 by an angle 7 to satisfy 
J22 = rn2. This rotation preserves J32 = m^. Moreover, because Ji^ + + </32 = 0, we get 
J22 = —nil — "^3 = "^2 automatically. 

To find the angle 7, we calculate the final vectors of applying the rotations -R(j3, 7) -R(y, /?) = 
i?(y,/3) i?(j3,7) to the vector in Eq. (|A2D . 



^ cos f3 cos 7 sin t]2 + sin /3 cos t]2 ^ 



Ji = Jl 



v 



sm 7 sm ri2 
sin (3 cos 7 sin ri2 + cos f3 cos 772 



J2 



^ — cos P COS 7 sin 771 + sin /3 cos 771 ^ 
— sin 7 sin rji 
sin /3 cos 7 sin rji + cos /3 cos rji 



(A5) 
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( sin/3 ^ 


cos /3 j 



Then impose the condition Ji^ = mi on the final vectors to determine 7, 



cos 7 



ji cos /3 cos ?72 — mi 



(A6) 



ji sin /3 sin 772 

In general, we find two solutions to Eq. (1A6|) . Let 7 represent the root of (]A6I) in the range 
[0,7r], and —7 the root in the range [— 7r,0]. 

We now lift the rotations -R(y, /3) -RQs, 7) up to an SU{2) rotation 



M(y,/3)M(z,7) 



(A7) 



(.Qg ^ /2 -e*T/2 sin /3/2 

and apply this spinor rotation flA7p to the reference spinors (\A3\i . This gives a set of spinors 



that project onto the final vectors in Eq. (]A5|) . This way, we have found one point p that 
satisfy the stationary phase conditions. The result is displayed in Eq. ( 134|) . 



Appendix B: The Calculation of the Hessian 

The stationary phase points consist of two 4-tori, so we pick eight transversal directions 
to calculate the Hessian. Let = r^^e*'^'''', s = 1, 2, 3, yU = 1, 2, be the complex coordinates 
written in polar coordinates. Introduce the new variables 

Zs = Zsl/Zs2 = RsC'"^ , (Bl) 

where Rs = rsi/r^2 and $^ = (psi - 4's2- 

We choose the eight directions transversal to the stationary phase points to be 9$^, 9$2, 
<9i?2) ^Ra^ ^ri2? dr22^ ^rz2- This is not an orthonormal coordinate system, so the change 
of coordinates generates a Jacobian factor J^^^ in the Hessian. 

Using the new variables Z^, the phase function / from Eq. (|24|) can be rewritten as 
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f = Yl [(-^'^ + 1^^^] + ^1 - ^3) + k2 ln(Z3 - 

s 

+k, HZ, - Z2) -Y^rl^Rl + 4j.lnr,2 - 



(B2) 



s2 ' 



Its first derivatives are displayed in Eqs. (]B3P - (IBlOp . Out of the 36 distinct second deriva- 
tives of /, there are 21 nonzero ones. These are hsted in in Eqs. fIBlip - fIBSip . 



5/ 
d^i 
dj_ 
d^i 
dl_ 
dRi 
dl_ 
dR2 

df_ 
dri2 
dj_ 

dr22 
dl_ 

drs2 



Zjj, — Zji 

-(j2 + "^2) + ^^^^^7 + 



1 
1 
1 



(jl + + ^^^l + 

{32 + ms) + + 

(j3 + "^3) + y ^^^y + 

— ^3 

4ji 







— Z2 




-kzZ2 ' 




Z\ — Z2 


1 


kzZx 




Z\ — Z2 


-k'iZ2 


— 9r^ R^ 

^' 22 ^2 


Z\ — Z2 


k2Z3 


— 2r32 -R3 


Zs-Zi 



-2(1 + R\)r^2 + — 
ri2 

-2(1 + i?2)r22 + — 
r22 

-2(1 + i?^)r32 + ^ 

?^32 



(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
(B8) 
(B9) 
(BIO) 



^2^1^3 
(^3 - ^1)^ 
/C3Z1Z2 

"(^1-^2)^ 



(9$1 9$2 

ay 

a$i dR2 

ay 

(9$2 5$2 



a$2 dRi (Zi - ^2)2 



fc3^1^2 



(z. 



L(^3-^l) 

(^1 - ^2)2 
2^2^! 

(^^^^ 

^1^2^3 
iZ2 - Z,)^ 

iA;3Z2e**i 



+ 



^3^26 



i<E>i 



(Zi 



+ 



fc3^1^2 
(^1 - ^2)2 



(Bll) 
(B12) 
(B13) 
(B14) 
(B15) 
(B16) 
(BIT) 
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d^2 dR2 



dR2dRs 

dRsdRs 

df 
dRi dri2 

df 
dR2 dr22 

df 
ORs dr32 

df 
dri2 dri2 

df 
dr22 dr22 

df 
drz2 dr32 



{Z2-Z3Y {Z^-Z2Y 
ik\Z2 



d^2 dRs 


(z 


2-^3)^ 




d'f 








dRidRi 




[Ri (Zs- 


-z,r 


d\f 


h 






dRidR2 




1 - ^2)2 




d'f 


k2 e^*i 




ORidRs 


(Z 


3 - Z,r 




d'f 




kiZ^ 




dR2dR2 




.-R2 {Z2 ■ 


-Zs)^ 


d'f 


ki e'*2 





+ ^ , „ — + 4 r 
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+ 



k^ e^*2 



(^2-^3)' 

k\ Z2 



+ 



k2Zx 



Rs{Z2-Zsr ' i?3(^3-^l)' ' ^ '\ 

= -AruRi 
= -4r22-R2 
= -4r32-R3 

= -2(l + i?? + ^) 
'12 

^-2{1 + Rl+^-^) 
'22 

= -2(l + i?3^ + ^) 

'32 



(B18) 
(B19) 
(B20) 
(B21) 
(B22) 
(B23) 
(B24) 
(B25) 
(B26) 
(B27) 
(B28) 
(B29) 
(B30) 
(B31) 



After simplifying the Hessian matrix through Gaussian ehminations on the last six rows 
and columns, which preserves the determinant, we find that the Hessian matrix has the form 
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H 





H12 



















\ 


H21 


H22 




























H33 




















































































Hm 




























H77 






























/ 



(B32) 



where 



-^11 ~ — ^"^^ 1 -^12 — H2I — T^zr — — — , H22 — TTT — ^TJ— , (B33) 

^33 = -K2 , ^44 = -4r^2 , ^55 = -4r^2 , (B34) 

//66 = -2(1 - i?? + ^) , //77 = -2(1 - i?^ + ^) , (B35) 

^12 ^22 

//88 = -2(1 -Rl+'^). (B36) 



At the stationary phase points, r^i = Vis + "^s' = Vis — "^s, s = 1, 2, 3, so the last three 
diagonal entries are constants. 



i/ss = i^66 = i^77 = -2(1 - i?^ + 2jVr22) 

-2 (^1 - -^^ + + 2j, 



(B37) 



-4. 



Taking the determinant, we find 



det /7 = 4^ 



^1 ^2 -Z^l -2^2 -Z^l 

(Z3 - z^y {Z2 ~ Z3) 

ki ks Z\ Z\ Z3 



+ 



^2 ^3 -Z^l -^2 Z3 



(Z3 - (Zi - Z2Y 



(B38) 



+ 



(Zi-Z2)2(Z2-Z3)2_ 

This determinant can be put into a more symmetrical form. 
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det H = 4:^ 



Zi Z2 



{Zi-Z2)iZ2-Zs){Zs-Zi] 
k2 h Zi {Z2 - Z-i) 



ki k2 Zs {Zi - Z2) 
(Z3 - Z^) {Z2 - Zs) 
ki k3Z2{Z3 - Zi) 



(^3 



-koZ- 



2^2 



Z3 — Zi 

k2Z3 



Zz — Z\ 

9[9i92 - gm] 



+ 



- ^1) {z, 

hZ2 
Z\ — Z2 
-kj.Z2 



[Zx — Z2) {Z2 — Z3 



kiZ-i 



Z2 ~ Z3 
-kiZi 



+ 



+ 



-k^Zi 
Z\ — Z2 
kzZx 



Z\ — Z2 ) \Z2 — Z3 Z\ — Z' 



(B39) 



where g corresponds to the factor in front, and g\-,g2igz-, gi correspond to the four factors in 
the four parenthesis, respectively. We now express them in terms of the variables z^^. We 
find 



9 = 4^7 ..^'^'^\. 7 (B40) 

{Z,-Z2)iZ2-Z3)iZ3-Z,) ^ > 

.6 ^11^12^21^22^31^32 (T}.AA\ 

Setting the first derivatives ( IB3p and flB4p to zero, and using jz^iP = jr + '^r, we can rewrite 
gi and (74 as 



— ^2^2 , ^3^2 /. , N ^2 _ 2:21 ^12 2;2i2;i22;ii f-r,Ac^\ 

gi = J ^ + J ^ = (Ji + mi)— = ziizx,—— = — , (B42) 

^3 — ^1 ^1 — ^2 ^1 Z22ZU Z22 

—kiZi k^Zi Zi _ Z\\Z22 Z11Z22Z21 /T-)^o\ 

£/4 = ^ ^ + ^ ^ = U2 + "^2)^ = 2:212:21 = • (B43) 

Z'2 — ^-3 Z-i — Zj2 ^2 Z\2Z2\ Z12 

Similarly, setting the first derivatives ( 1B3I) and ( ]B4I) to zero, and using 2ji = k2 + k^, 
2^2 = ki + ks, and |zr2p = jr — ^r, rewrite g2 and g^ as 



^1^3 , —^3^1 kiZ2 —k3Z2 , 1 \ 

g2 = ^ ^ + ^ ^ = ^ ^ + ^ — ^ - (fci + K3) 

^2 ^ ^3 ^1 ~ ^2 ^2 ^ ^3 ^1 ~ ^2 

= (j2 + "^2) - 2j2 = -(j2 - m2) 

= -Z22Z22 , (B44) 

^2-^3 , —k2,Z2 ^2^1 —^3^1 , /, , 1 \ 

gs = ^ ^ + ^ ^ = ^ ^ + ^ ^ + [h + ks) 

^3 ~ ^1 -^1 ^ ^2 ^3 ~ -^1 ^1 ~ ^2 

= -(ji + + 2ji = (ji - nil) 

= Z\2Z\2 ■ 
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Using the values of gi, g2, g^, g4, and g above, we find 



detif = g{gig2 - g^Qi) = 5' (-^21 ^12 2:11 ^22 + 2;ii2;22^i2^2i) 

= -AtgA,, (B45) 
where we have used Eq. fl30|) - fl32l) in the third equahty. In the last equality, we have defined 

A. = l{JlxJ2y - J2xJly) ■ (B46) 

It is the area of the triangle in the angular momentum space, projected onto the x-y plane. 
Finally, the Hessian is equal to det H times the Jacobian Ylsfj. ^"sf^- The result is 



Hessian 



2;ii2;i22;2i2;222;3i^32 



(2:11^22 ~ Zi2Z2l){z3iZi2 — ^32^11) (^21^32 ~ 2:222:31) 

(-^)4^A. ^ 

ZiiZi2Z2lZ22Z3lZ32{ZiiZ22 — Z12Z21) {Z31Z12 — Z32ZU) {Z21Z32 — Z22Z31) 



. (B47) 



Note that the factor —i is responsible for the Maslov index 1/ 



I = e 
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